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Abstract 

We derive the topological susceptibility to the one-loop order in the chiral effective theory of 
QCD, for an arbitrary number of flavors. 
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I. INTRODUCTION 



In Quantum Chromodynamics (QCD), the topological susceptibility (xt) is the most 
crucial quantity to measure the topological charge fluctuation of the QCD vacuum, which 
plays an important role in breaking the Ua_{^) symmetry. Theoretically, Xt is deflned as 



Xt = I d'x {p{x)p{0)) , (1) 

where 



p{x) = ^^e^^xatT[F^Ax)Fxa{x)], (2) 



is the topological charge density expressed in term of the matrix- valued fleld tensor F^iy. 
With mild assumptions, Witten l|] and Veneziano 2(] obtained a relationship between the 
topological susceptibility in the quenched approximation and the mass of rj' meson (flavor 
singlet) in unquenched QCD with Nj degenerate flavors, namely. 



(quenched) 



2Nf 



where Ft^ ~ 93 MeV, the decay constant of pion. This implies that the mass of rj' is 
essentially due to the axial anomaly relating to non-trivial topological charge fluctuations, 
which can turn out to be nonzero even in the chiral limit, unlike those of the (non-singlet) 
approximate Goldstone bosons. 

Using the Chiral Perturbation Theory (ChPT), Leutwyler and Smilga 0, 4] obtained the 
following relations in the chiral limit 

Xt=^(— + —y\ {Nf = 2), (3) 



rriu TTid. 
1 1 1^-1 



Xt=S — + — + — , {Nf = 3), (4) 

where m^, m^, and rris are the quark masses, and S is the chiral condensate. This implies 
that in the chiral limit (m^ —>■ 0) the topological susceptibility is suppressed due to internal 
quark loops. Most importantly, ([3]) and (jl]) provide a viable way to extract S from xt in the 
chiral limit. 

From ([T]), one obtains 

Xt=^-^, Qt^J d'xpix), (5) 



where Q is the volume of the system, and Qt is the topological charge (which is an integer 
for QCD). Thus, one can determine Xt by counting the number of gauge configurations for 
each topological sector. Furthermore, we can also obtain the second normalized cumulant 

c. = ~[{Qt)-mr], (6) 

which is related to the leading anomalous contribution to the rj' — rj' scattering amplitude 



in QCD, as well as the dependence o: 



the vacuum energy on the vacuum angle 6. (For a 



recent review, see for example, Ref. js] and references therein.) 

However, for lattice QCD, it is difficult to extract p{x) and Qt unambiguously from the 
gauge link variables, due to their rather strong fluctuations. 

To circumvent this difficulty, one may consider the Atiyah-Singer index theorem 6| 

Qt = — n_ = index (D), (7) 

where n± is the number of zero modes of the massless Dirac operator V = 7^(9^ + igA^) 
with ± chirality. Since V is anti-Hermitian and chirally symmetric, its nonzero eigenmodes 
must come in complex conjugate pairs with zero chirality. Thus one can obtain the identity 



n I — n_ 



m J d'^x tr[y^{T) + m) ^(x,x)], (8) 



by spectral decomposition, where the nonzero modes drop out due to zero chirality. In view 
of ([7]) and ([H]), one can regard pt{x) = mgtT['y^{V + mg)~^{x, x)] as topological charge density, 
to replace p{x) in the measurement of Xt- 

Recently, the topological susceptibility and the second normalized cumulant have been 
measured in unquenched lattice QCD with exact chiral symmetry, for Nf = 2 and Nf = 2 + 1 
lattice QCD with overlap fermion in a fixed topology [3, Iq], and Nf = 2 + 1 lattice QCD with 
domain- wall fermion [3] . The results of topological susceptibility turn out in good agreement 
with the Leutwyler-Smilga relation, with the values of the chiral condensate as follows. 

E^(2 GeV) = [245(5)(12) MeV]^ (A^/ = 2), Ref. [7], 
E^(2 GeV) = [253(4) (6) MeV]^ (A^^ = 2 + 1), Ref. [8], 
E^(2 GeV) = [259(6) (9) MeV]^ {Nf = 2 + 1), Ref. [9]. 

These results assure that lattice QCD with exact chiral symmetry is the proper framework 
to tackle the strong interaction physics with topologically non-trivial vacuum fluctuations. 



Obviously, the next task for unquenched lattice QCD with exact chiral symmetry is to 
determine the second normalized cumulant C4 to a good precision, and to address the question 
how the vacuum energy depends on the vacuum angle 9 and related problems. Theoretically, 
it is interesting to obtain an analytic expression of C4 in ChPT, as well as to extend the 
Leutwyler-Smilga relation to the one-loop order. In this paper, we derive the topological 
susceptibility to the one-loop order in ChPT, for an arbitrary number of flavors. 

The outline of this paper is as follows. In Section 2, we review the derivation of topological 
susceptibility Xt at the tree level of ChPT, and also derive the second normalized cumulant 
C4 at the tree level, and discuss its implications. In Section 3, we derive Xt up to the one- 
loop order in ChPT for an arbitrary number of flavors. In Section 4, we conclude with some 
remarks, and also present the case of 2 + 1 flavors, in which only the one-loop corrections due 
to the u and d quarks are incorporated. In the Appendix, we present a heuristic derivation of 
the counterpart of the Leutwyler-Smilga relation in lattice QCD with exact chiral symmetry. 



II. TOPOLOGICAL SUSCEPTIBILITY AT THE TREE LEVEL OF CHPT 

Before we proceed to derive xt to the one-loop order in ChPT, it is instructive for us to 
recap the derivation of Xt at the tree level 0, 4|. 

The leading terms of the effective chiral lagrangian for QCD with Nf flavor at = [l^ 
are the kinetic term and the symmetry breaking term, 

£(2) = + 4i = ^Tr(a^f/a'^f/t) + |Tr(A^t/t + UM^), (9) 

where U{x) = exp{2i(l)"'{x)t"- / F.„} is a group element of SU{Nf), M. is the quark mass 
matrix, is the pion decay constant, and S = {ipip)vac is the chiral condensate of the QCD 
vacuum. 

On the other hand, the partition function of QCD in the 6 vacuum can be written as 

ZN,{e) = Y.e-''^'ZQ, (10) 

Q 



where 



5 1 _L^,5 



Zq = J [ciA^]e-^«[^-l det {j,D^ + -^M + -^^M^j 

= /MA,]e-«'-^lndet(A^ + A^t^).|'^'^^'^' ^ ^ °' 
^ k det{M)-'^, g < 0, 



where Sg is the action of the gauge field, and A^'s are non-zero eigenvalues of the massless 
Dirac operator in the gauge background. Thus the physical vacuum angle on which all 
physical quantities depend is 6'phys = + argdet(A^) rather than 9. Also, the 6'-dependence 
of Z^^iO) always enters through the combinations Ade^^^'^f and A^^e"*^/^^. It follows that 
for ^ 7^ 0, the symmetry breaking term in the chiral effective lagrangian can be written as 



^S. = S Re 



Defining the vacuum energy density 



1 



(11) 



(12) 



then the topological susceptibility Xt © and the second normalized culmulant C4 ([6]) can be 
expressed as 

<9^evac(A^,6 



Xt 



(13) 
(14) 



For small quark masses {L <^ ^^^), the unitary matrix U does not depend on x^. Thus 
only the symmetry-breaking term survives in ([9]), and the partition function becomes 

Z^^^e) = J dUexp {q S Re [TiiMe'^^^fU^)] } , (15) 

where Q = L^T is the space-time volume. Without loss of generality, the unitary matrix U 
can be taken to be diagonal 



t/ = diag (e^°S e^"% 



0, 



where the last constraint follows from the unitarity of U. Also, we can choose the mass 
matrix to be diagonal Ai = diag(mi, . . . , m^v^). Then we have 



Re 



^31 



where (pj = 9/Nf — aj, and J2j (pj = 0. 

Now, we consider a sufficiently large volume Q satisfying rrijllQ ^ 1, then the group 
integral in the partition function f[T^ is largely due to the U which minimizes the minus 
exponent of the integrand, i.e.. 



min {-Re [TiiMe'^/^fU^)] } = min <^ - ^ rrij cos 



e. 



(16) 



i=i 



For Nf = 2, this amounts to minimize the function 

—mi cos(0i) — 1712 cos{9 — (pi), 



where the constraint 0i + 02 = ^ has been used. A simple calculation gives the minimum, 



— Y mf + 1712 + 2mim2 cos 9. 



Thus the partition function is 



ZNf{9) = Zoexp jfiS^mf + + 2mim2 cos 6'| , mEi7 ^ 1, 



which gives the vacuum energy 



evac{9) = €o - T^^ml + ml + 2m„mrfCos6', 

where eo is the additive normalization constant corresponding the normalization factor Zq 
in the partition function. From f|T71) . we obtain the topological susceptibility 



Xt 



09^ 



'rriu + ma 



(17) 



Furthermore, the second normalized cumulant is 



C4 



89^ 



m^ma 



3S- 



mlml 



\ mt 



rriu + ma 



(18) 



=0 mu + ma {m^ + m^Y 

which has not been discussed explicitly in the literature. The vital observation is that the 
ratio of Xt and C4 is 



C4 
Xt 



-1 + 



(19) 



(m„ + mdY 

which goes to —1/4 in the isospin limit m„ = m^. This seems to rule out the dilute instanton 
gas/liquid model [ll, Q, Q which predicts that C4,l_Xt = — 1. Moreover, recent numerical 
results of C4,/xt from quenched lattice QCD [isl, 16] and unquenched lattice QCD 0, 
are consistent with the prediction of ChPT. 

Next we turn to the case Nf > 2. Then there is no analytic solution to the minimization 
problem (fT6|) . However, for the purpose of obtaining the topological susceptibility, one may 
consider the limit of small 9 (and 0j's) because U = 1 gives the minimal vacuum energy at 
9 = 0. Since Xt only depends on the curvature of evac(^) around = 0, this approximation 



6 



would give the exact result of Xt (at the tree- level). To the order of 6'^, the minimization 
problem f|T6|) becomes 



Nf 



mm 



COS( 



mm 



I Nf 



Nf 

1=1 



9. 



Now introducing the Lagrange multiplier A to incorporate the constraint J2i 0i = then 
the minimization problem amounts to solving the equation 

'Nf 



d 



1^ ,o 1 



-^m,-02 - - A I E<^J - ^ 



mi(f)i - -mj0,^ - A = 0. 

D 



Setting (f)i = + {m') (where ai and 03 are parameters), and using J2i 0i = ^) we can 



solve for oi and 03, and (pi to the order of 9^, 



m ^ 9^ 
m,- D 



m 
rrii 



m 

TTli 



Nf / - \ 3' 

' m \ 



E 



+ o(e^). 



where m = (j^fji is the "reduced mass" of the Nf quark flavors. Keeping the 

exponent of the partition function to the order of we have 



Nf 



ZA./0) = Zoexp|-r]E 
and the vacuum energy density is 



i=l 



u=i 



24 



+ 0{9^) 



(20) 



Nf 



-1 



02 



iV 



e...(0) = eo + S|5:i-) y-SE"^^- 



0^ 



i=l 



E-l ^ + ^(^^ 



It follows that the topological susceptibility is 



and 



Xt 



Nf ^ 



-1 



(21) 



C4 



E 

i=l 



'Nf . 

E — 



(22) 



which generalize Eqs. f[T7|) and flTSl) to an arbitrary number of flavors. In particular, for 



and 



(I 1 1 
C4 = -S -^ + -3 + -^ — + — + — • (24) 



Nevertheless, these two formulas seem unnatural, since the strange quark is much heavier 
than the up and down quarks. Thus a plausible chiral limit is to take mu,d 0, while keeping 
rris fixed. Consequently, the condensate of the strange quark (ss) must be different from E, 
and it should also enter this formula. In the Appendix, we present a heuristic derivation of 
the counterpart of (1231) in lattice QCD with exact chiral symmetry, which takes into account 
of the difference between (ss) and E, as given in Eq. (l52l) . 



III. TOPOLOGICAL SUSCEPTIBILITY TO THE ONE-LOOP ORDER OF CHPT 



To the one-loop order of ChPT, one has to include C^'^'' |lO| at the tree level as well as the 
one-loop contributions of C^'^\ In 1984, Gasser and Leutwyler {lo| considered the low-energy 
expansion, where both p and A4 are assumed to be small but M./p^ can have a finite value, 
such that the value of M^/p"^ can be fixed. In this case, the external sources a^(a;) and p{x) 
can be counted as order of $, and f^(x) and s{x) — Ai as order of Gasser and Leutwyler 
showed that at the one-loop order, the chiral effective action can be written as 

W = Wt + Wu + WA + 0{<i>^), (25) 

where Wt denotes the sum of tree diagrams and tadpole contributions (of order $^), Wu the 
unitarity correction (of order $^), and Wa the anomaly contribution (of order $^). Because 
the 6 dependence enters the Lagrangian only through Ad, we can count Xt as order of 
thus for the evaluation of topological susceptibility to the one-loop order, and it suffices to 
consider Wt only. 

□ 

Moreover, Gasser and Leutwyler showed that the pole terms due to the one-loop 
contributions of C^'^^ can be absorbed by the low-energy coupling constants of C^^\ and Wt 
is given by 

Mi 



8 



where M|,'s are the squared meson masses, app corresponds to the kinetic term which can 
be dropped in the hmit of small quark masses, app corresponds to the symmetry breaking 
term, 



^pp 



lTr({Ap,Al,}(xtf/ + t/tx))-M| 
and is just C^^"* with renormalized low-energy coupling constants. 



(27) 



Ci^)=L\[Ti[D^U{D''U)' 



D^UiDMV Tr D^'UiD'^Uy 



+LITT 
+LITT 



D^U{D^'UyDM{D''Uy 
D^U{D^Uy] Tr (xU^ + Ux 



+LITT iD^UiD^UYixU^ + Ux^)] + [Tr (xU^ + Ux 



+U, [Tr [xU^ - Ux^)] + L^Tr iUx^Ux' + xU^U^ 



(28) 



Here x = 2(S/F^)7V/i = 2BqM.^ Ap's are the generators of SU{N) in the physical basis, 
{U'{fisub),i = I5 ■ ■ ■ ) 10} are renormalized low-energy couphng constants, and the last two 
contact terms (with couplings HKfisub) and Hl{iisuh)) are the counter terms required for 
renormalization of the one-loop diagrams. 

For small quark masses (L ^ "^^^), the unitary matrix U does not depend on x^, thus 
the term involving app in (!26|) can be dropped. 

Next we consider the term with a^p in fl26l) . Using the formula 

r 



A{N} 



Nf 



I, 



p 

we obtain its contribution to the chiral effective lagrangian. 



(29) 



EReTr(Xt/f) 



NfF 



2{N} 



EReTr({A.,At}^^t)^ln^ + -|^ln 



Ml 



Mi 



TT P 



l^sub 



Ml 

f^sub 



(30) 



For small quark masses (L <^ ^tt^)? ^^e unitary matrix U does not depend on x^, so only 
the sixth, seventh, and eighth terms in are relevant to the partition function. For 9^0, 



9 



6 enters the chiral effective lagrangian through the combinations Aie^^^^f and Ai'^e~^^^^f . 
Thus these three potential terms can be written as 

2 r , -o/AT j-,n2 



4SoReTr(A^e*^/^-ff/"^) + iABolmTi{Me'^/^fU^) + SL^E^ReTr (Me'^^/^fU^f 



Without loss of generality, we can take U and Ai to be diagonal, 

Nf 

U = diag(e^"^e^°^...,e^°-/), E = 0> 



(31) 



M = diag(mi, . . . , mATj,), 



therefore the contributions of fl30l) and fl3Tl) become (dropping the terms without ?7 depen- 
dence) 



Nf 



E X] m,- cos E E l^^^' ^p},-, "^i cos cj)j^ In ^ 



Nr 



Nf 



Nf 



+ IQBqLq E "^i cos ( 



165^^^ I E sin (pj + SB^Ll ^ cos : 



where 0j = ^/A^j — aj, and X]j 0j = ^• 

Again, we use small 6 (small 0/s) approximation and keep terms up to the order of 
then the evaluation of the integral in the partition function in the limit mjf2S ^ 1 amounts 
to minimizing the generating functional 



mm 



Nf Nf 



' Nf 



Nf 



_ _ (32) 

i=i j=i 

with the constraint J2j 'Pj = ^- We introduce the Lagrange multiplier A to incorporate the 
constraint in finding the minimum. For simplicity, we define 

Nf 



Ai = —rrii — 



E{ap,a1.} 



3 J 167r2 /if^^ 



^Mn ^ + 16i?o' I Llm, Y.m. + L.. 



TT P 



1=1 



Bj = ABoir^y/^nij. 



Then the minimization problem amounts to solving the equation 



d 



Nf 



'Nf 



'Nf 



E^.0|+ E^. 



10 



which gives 



Nf 



Ai(j)i + BiiYl BjC^j 
Defining {T)ij = 2Ai5ij + 2BiBj, ([33]) becomes 



A 
2' 



(33) 



$:(T),,0,- = A, z = l,...,iV;. 



(34) 



Thus we can obtain A using the constraint 

Nf Nf Nf 

j=i j=i i=i 



X = 9 



'Nf Nf 

j=i 1=1 



(35) 



Now multiplying Eq. fl33|) with 0j and summing over i, we obtain 

21 



mm 



TV, 



Nf 



i=i Vi=i 



A^ 

Y' 



(36) 



which can be used to simplify (152]) to 

xe _ 9^ 



Nf 



E (T-^).. 



(37) 



where the last equality follows from (135]) . Finally, the partition function of QCD with Nf 
flavors to the one-loop order of ChPT in the limit mSi7 ^ 1 is equal to 

-1' 



Zn^ (9) = Zq exp < 
and the vacuum energy density is 

evac(6') = Co + 



.9' 



Nf 

E (T-^: 



(38) 



9^ 



■ Nf 

E (T-^; 



(39) 



Thus the topological susceptibility to the one-loop order of ChPT is 

-1 



Xt 



89^ 





" Nf 




E (T-^)., 


61=0 





(40) 



To simplify the expression of topological susceptibility, we rewrite the matrix T as 

(T),, = 2AAj + 25^5,- = S(M + T')ij, (41) 
11 



where 



Ml 



J 777. J , 



fc=l 



and 



32 



sub J 



(42) 



(43) 



Since the eigenvalues of the real and symmetric matrix A4 ^/^T'A^ are much less 
than one in the chiral limit, we can use the Taylor expansion 

and obtain 



Xt 



■ Nf 



^/ fT'V- 



4^2 



m\ Ml ^ Ml 
' m 



su6 



+Ks Y.mi + NfiNfKr + Ks)m 



(44) 



i=l 



where m = (YI^Ji m^^^ , and all terms proportional to Kf or KiKj have been dropped. 
Equation (jH]) is the main result of this paper. 

For Nf = 2, there are three mesons, tt"*", 7r°, and tt^. If we take their masses to be the 
same and use (l29l) . we obtain 



Xt 







-1 r 




ma J 





3 

^ln-^ + i^6(m„ + mrf) 



2F2 167r2 ^2^, 



(45) 



Next we turn to the case Nf = 3. Taking the eight pseudoscalar mesons with non 



degenerate masses, we obtain 

1 



Xt 



Sm< 1 



2F2 



m ^ BQlrrii + rrij) ^ BQ{mi + rrij) 



m,; rrij ^ 



m 



m \ M% , M\ I f m 



+ — + — -7^1n^ + - — + — +4— 



m 



//2 

H'sub 



m 



(46) 



12 



where m = (m^^ + ^ + m^^^ , and Bq = 

IV. CONCLUDING REMARK 

In this paper, we have derived the topological susceptibility to the one-loop order in 
ChPT, in the limit mEVi > 1, for Nf = 2 [Eq. (US])], Nf = 3 [Eq. (^], and an arbitrary 
number of flavors Nf [Eq. (jSj)] respectively. 

For Nf = 3, since the mass of the strange quark is much heavier than the masses of u 
and d quarks, it seems reasonable just to incorporate the one-loop corrections due to the u 
and d quarks. Then, for Nf = 2 + 1 {u and d quarks to the one-loop order, and s quark at 
the tree level), the topological susceptibility becomes 



xt = ^U — + —) 



?, Ml ^ Ml , 
2F2 167r2 ^^^^ 



+ 1^1 ■ (47) 



This supplements (l46l) for the case A^/ = 2 + 1. 

Now the trend of unquenched lattice QCD simulations is to include the charm quark. 
Thus it is also interesting to include the case N f = 2 + 1 + 1, with both s and c quarks being 
kept at the tree level. Then the topological susceptibility is 



Xt = ^{{ — + — 



3 Mf . M2 



^ In -^^ - K^^rriu + m^) 



2F2 levr^ 

2i2K, + Ks) "^""^^ 



+ — + -I (48) 



rriu + md. 

In view of the one-loop results of Xt, [Eqs. ( l45l) . ( l46|l . ( l47l) . and ( l48l) ]. it would be 
interesting to see whether the Xt measured in lattice QCD with exact chiral symmetry 
would agree with the prediction of ChPT. Most importantly, these one-loop formulas provide 
a viable way to determine the low-energy constants F^, Lq^L-? and Lg, in addition to the 
chiral condensate S which has already been determined {3, S, Is] using the formula of Xt at 
the tree level fl23|) . 

Finally, we turn to the second normalized cumulant C4. At this moment, we only have a 
formula of C4 (122|) at the tree level. For Nf = 2, the ratio c^/xt = —1/4 in the isospin limit 
{rriu = md) seems to rule out the instanton gas/liquid model which predicts that c^/xt = — 1- 



13 



Obviously, it would be interesting to derive a formula of C4 to the next (non- vanishing) order 
in ChPT. 



Appendix 

In this Appendix, we present a heuristic derivation of the relationship between topological 
susceptibility, chiral condensate, and the quark masses, for an arbitrary number of (non- 
degenerate) flavors, in the framework of lattice QCD with exact chiral symmetry. Our 
derivation generalizes that presented in Ref. l3] with degenerate flavors. 

Consider the flavor- singlet pseudoscalar v( 

f 1=1 

Its correlator at zero momentum is 



1 



Nf 



H H {.<ii{x)ibqi{x)qj{y)imj{y)) 



[dU]deiD{m)e 



X 



Nf 



Nf 



i=l 



vlzn 



\ r 1 
—2 / [dU] det D(m)e-^«[^1 ^ — Tr(/^c + 



nii 



Nf 



n< — n_ 



i=i 



(49) 



where Sg[U] is the gauge action, 

detL'(m) = \{det[{D^ + mi){l + rD^)-\ 



Nf 



i=l 



[dU] det D{m)e 



-Sgiu] 



and the identity 



Tt[{D, + m)-SPc + m)-S5] = -Tr(D, + m)~\ 



m 



has been used in the last equality of (H9|) . Here Z)c = D{1 — rD) ^ is the chirally symmetric 
Dirac operator of a Ginsparg- Wilson Dirac operator D satisfying D75 -|- 75D = 2rD'-)^D. 
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Now taking the thermodynamic hmit (f2 oo), and then the chiral hmit (m^ 0), (H9 
gives 



1 /'^ 1 \ f„ 1 



G.'(0) = 1^|E-MS-|E-|X*^ (50) 

where 



E = hm hm ^{Tr{Dc + rrii)-^) , 
^* = hm hm ^{{n+ -n^f). 



If 1]' stays massive, then its propagator G^/ oc m^,^ must be non-singular. This imphes that 
the coefficient of the singular factor (j^fji ^) in (150|1 behaves like C(m), i.e., 



1 

1 



i=i 

which agrees with the Leutwyler-Smilga relation. 
For Nf = 2 + 1 with fixed m^, (!50|) is modified to 

where 

{ss) = lim ^(Tr(£)c + 

In the limit m„^rf 0, the coefficient of the singular factor 2/mu,d in (I^Tl) behaves like 
0{mu,d), i.e., 

/ E E / 1 1 1 \ 

which provides a more physical relationship (between xt, S, (ss), and the quark masses) 
than the Leutwyler-Smilga relation (jlj), since it reduces to (jll) only in the (unphysical) limit 
(ss) = E. 

Now it is straightforward to generalize fl52|) to the case Nf = 2 + 1 + 1, 

/ E E (ss) (cc)\ / I 1 1 1 
rud rus j Vm„ 

where 

(cc) = lim \-l(lx(Dc^m^'^). 

It would be interesting to determine E, (ss), and (cc) with the data of in lattice QCD 
with exact chiral symmetry. 
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